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NOISE TRANSMISSION BY VISCOELASTIC SANDWICH PANELS 

Rimas Vaicaitis* 

Langley Research Center 


SUMMARY 

This report presents an analytical study on low frequency noise transmis- 
sion into rectangular enclosures by viscoelastic sandwich panels. The dimen- 
sions of these panels are typical of aircraft fuselage skin constructions. The 
report demonstrates that these panels can effect significant noise reduction. 

The analysis considers two limiting cases of the core stiffness: (1) soft com- 

pressible cores with dilatational modes included and (2) hard incompressible 
cores with dilatational modes neglected. Sandwich panels with soft viscoelastic 
cores exhibit noise transmission characteristics similar to those of double wall 
elastic panels except in the frequency range where dilatational modes occur (300 
to 600 Hz). At these frequencies noise can be reduced by as much as 30 dB more 
with viscoelastic panels. Sandwich constructions with hard cores show signifi- 
cant noise reduction advantage over the elastic panels for the whole frequency 
range considered in this study (0 to 1000 Hz). Numerical results contained in 
this report include response and noise transmission characteristics of elastic 
and viscoelastic panels for several geometric and material parameter 
configurations . 


INTRODUCTION 

The information available in the literature and from ongoing research pro- 
grams on interior aircraft noise indicates that noise in many aircraft exceeds 
acceptable comfort limits. Propeller driven aircraft where maximum noise inten- 
sity occurs at low frequencies are especially blameworthy. Since acoustic 
absorption materials used in aircraft constructions are not very effective in 
reducing interior noise at low frequencies (ref. 1), new means of providing 
noise attenuation at low frequencies need to be established. Interior aircraft 
noise in the structural resonance range is strongly controlled by the vibra- 
tional characteristics of the fuselage skin panels. Past studies have demon- 
strated that a viscoelastic material sandwiched between two elastic plates is a 
very efficient way of dissipating vibrational energy (refs. 2 to 9). Replacing 
some of the elastic skin panels with viscoelastic sandwich constructions should 
achieve significant amounts of noise reduction. Since available information on 
noise transmission characteristics of viscoelastic sandwich panels which are 
suitable for aircraft constructions is very limited, an analytical study of this 
subject is undertaken in this paper. The actual aircraft fuselage construction 
is too complicated for a detailed mathematical treatment; therefore, a simpli- 
fied analytical model has been constructed. The interior acoustic enclosure 
into which noise was transmitted is taken as a rectangular box. One wall is 
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flexible, and the remaining walls are rigid. (See fig. 1.) The dimensions of 
the flexible panel were chosen from typical aircraft skin panels, and the depth 
of the acoustic enclosure corresponds to about one-half of the aircraft cabin 
width. 

In this paper the noise transmission into the enclosure is analyzed by 
solving the linearized acoustic wave equation for the interior noise field and 
the plate vibration equation for the viscoelastic sandwich panel vibrations. 

The acoustic equation is coupled to panel vibrations through the time-dependent 
boundary conditions. The solution to this system of equations is obtained by 
using modal expansions and a Galerkin-like procedure. Since the boundary condi- 
tions for these equations are time dependent, the commonly used method of sepa- 
ration of variables cannot be applied to this system (ref. 10). The time depen- 
dence, however, is removed by splitting the solution into two parts: a solution 

corresponding to a nonhomogeneous differential equation with homogeneous bound- 
ary conditions and a solution on the boundary. Following this procedure, a 
Fourier series solution is developed which converges rapidly not only in the 
interior acoustic space but also on the boundary. The solution possesses con- 
tinuous derivatives up to the second order and satisfies the given partial dif- 
ferential equation and all the boundary conditions. It can be demonstrated that 
such a series solution is uniformly convergent (ref. 11). These series have a 
computational advantage over the nonuniformly convergent series which usually 
converge very slowly. 

The governing differential equations for the vibration of the viscoelastic 
sandwich panels due to the prescribed external pressure are developed for two 
limiting cases of core stiffness. In both of these cases the core material is 
assumed to be isotropic and the elastic plates are perfectly bonded to the core 
so that no displacement discontinuities are present. In the first case, the 
viscoelastic material is taken to be very soft so that bending and shearing 
stresses can be neglected, and the core acts merely as a viscoelastic spring. 

For the sandwich construction with this soft core, the flexural vibration modes 
are governed by the stiffness of the two face plates, and the out-of-phase dila- 
tational modes are controlled by the stiffness characteristics of the viscoelas- 
tic spring. In the second case, a stiff viscoelastic core material is used 
where the bending and shearing strains in the core are important. Because of 
the stiff character of this core material, dilatational modes occur at high fre- 
quencies. By limiting the analysis to low frequencies (below 1000 Hz), the core 
material can be assumed to be incompressible in this frequency range, and the 
out-of-phase motions of the two elastic face panels associated with the dilata- 
tional modes can be neglected. Furthermore, the stiff core model is simplified 
to satisfy the following assumptions: the shearing strains across the depth of 

the face plates are small, the stresses in the core material parallel to the 
plate surface are negligible, and the in-plane "rotary" inertia effects are not 
included. Since the analysis is limited to low frequencies and relatively thin 
viscoelastic cores, these assumptions are believed to be valid. For the stiff 
core panel, the governing equations of motion and the solution for panel vibra- 
tions are obtained from reference 2. 

This report contains numerical results for a rectangular enclosure with 
acoustically hard walls. It is assumed that one face is covered with a flexible 
viscoelastic sandwich panel as shown in figure 1 ; otherwise , the enclosure is 
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rigid. Results include panel response spectral densities, noise transmission 
characteristics, and overall interior noise levels. A direct comparison of the 
results between elastic and viscoelastic sandwich panels is shown. 

SYMBOLS 

The units used for physical quantities defined in this paper are given in 
the International System (SI) of Units. Correlations between this system of 
units and U.S. Customary Units are given in reference 12. 

A conversion factor 0.000645, see reference 12 

A^ j variables defined in equation (11) 

nf 

a plate length, m 

a-|,a 2 variables defined in equation (30) 

Bp flexible boundary 

B ijk pressure modal coefficients, N/m 2 

Bp rigid boundary 

b plate width, m 

0^ displacement modal coefficients of top face plate, soft core visco- 

elastic sandwich panel, m 

(Cp)ij pressure coefficients, N/m 2 

c speed of sound in cavity, m/sec 

ci damping coefficient of top elastic plate, N-sec/m3 

C 2 damping coefficient of bottom elastic plate, N-sec/m3 

D = Eh3/I2(l - v 2 ), plate stiffness, N-m 

Dmn displacement modal coefficients of bottom face plate, soft core 

viscoelastic sandwich panel, m 

D*] = E -j in-j 3/ 12(1 - V -j 2 ) , top face plate stiffness, N-m 

D 2 = E2h23/12(1 - V 2 2 ), bottom face plate stiffness, N-m 

d cavity depth, m 

E modulus of elasticity of plate, N/m 2 

E-| modulus of elasticity of top face plate, N/m 2 
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e 3 

e ij 

e k 

F ijk 

f 

f mn 

G 

G ij 

GO 

H ijk 

H mn 


h 

h 1 ,h2> h 3 


modulus of elasticity of bottom face plate, N/m 2 
modulus of elasticity of viscoelastic core, N/m 2 
variables defined in equation (7) 
variables defined in equation (19) 
variables defined in equation (18) 
frequency, Hz 

panel modal frequencies, Hz 

complex shear modulus defined in equation (36), N/m 2 
variables defined in equation (6) 
shear modulus of viscoelastic core, N/m 2 
variables defined in equation (54) 

frequency response function of sandwich plate with hard viscoelastic 
core, m3/N 

frequency response function of bottom face plate in sandwich construc- 
tion with soft core, m3/N 

frequency response function of top face plate in sandwich construction 
with soft core, m^/N 

plate thickness, m 

thickness of top face plate, bottom face plate, and viscoelastic core, 
respectively, m 



indices 


i imaginary unit (-1)^ 2 

K mn variables defined in equation (40) 

kg spring stiffness of soft viscoelastic core, N/m3 

L ijmn variables defined in equation (55) 

NT noise transmission, dB; see equation (59) 

OASPL overall sound pressure level, dB 
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p sound pressure in cavity, N/m 2 

PO reference pressure, 0.00002 N/m 2 

p r external random input pressure, N/m 2 

Qmn generalized random forces, see equation (28) 

R x spatial correlation coefficient corresponding to x-coordinate 

Ry spatial correlation coefficient corresponding to y-coordinate 

r*i variables defined in equation (44) 

Smnrs cross-spectral densities of generalized random forces, (N/m 2 ) 2 /Hz 
Sp spectral density of sound pressure in interior, (N/m 2 ) 2 /Hz 

Sp r spectral density of input pressure, (N/m 2 ) 2 /Hz 

S r cross-spectral density of input pressure, (N/m 2 ) 2 /Hz 

S w panel deflection response spectral density, m 2 /Hz 

t time, sec 

Umn>^mn (Wjjjn modal displacement coefficients corresponding to x-, y-, and 

z-coordinates, respectively, m 

u,v,w displacement components corresponding to x-, y-, and z-coordinates, 
respectively, m 

Wi transverse displacement of top face plate, m 

»2 transverse displacement of bottom face plate, m 

^ijk acoustic modes 

x,y,z spatial coordinates, m 

Z d variables defined in equation (52) 

variables defined in equation (13) 
a ijk acoustic cavity modal damping coefficients 

plate modal damping coefficients 
B damping loss factor in viscoelastic core 

Ynjn variables defined in equation (29) 
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^mn variables defined in equation ( 4 1 ) 

®ijk variables defined in equations (53a) and (53b) 

^mn variables defined in equation (29) 

y variable defined in equation (35) 

v Poisson’s ratio 

v 1> v 2> v 3 Poisson's ratio of top face plate, bottom face plate, 
and viscoelastic core, respectively 

p material density of elastic plate, kg/m3 

p a air density in enclosure, kg/m3 

Pl,p 2 ,P 3 material densities of top face plate, bottom face plate, 
and viscoelastic core, respectively, kg/m3 

a 1,02 variables defined in equation (30) 

t variable defined in equation (43) 

$ij variables defined in equation (14) 

<j) mn structural modes of panel 

4>mn variables defined in equation (42) 

a) frequency, rad/sec 

w ijk cavity modal frequencies, rad/sec 

structural modal frequencies, rad/sec 

Superscripts: 

b bottom 

t top 

Fourier transform according to equations (20a) and (20b) 
* complex conjugate 

Dots over symbols denote time derivatives. 
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INTERIOR ACOUSTIC PRESSURE 


The rectangular enclosure shown in figure 1 occupies a designated volume 
V = abd. The wall at z = 0 is flexible whereas the remaining walls are acous- 
tically rigid. The pressure p inside the enclosure is determined from the 
linear acoustic wave equation 

V 2 p = — p ( 1 ) 

c^ 

where V2 = 32/3 x 2 + 32/3y2 + 32/3 z 2 anc j the boundary conditions to be satis- 
fied are 


9p 

3n 


0 


on Bp and 


9 P 

— s 

3n 


-p a w 


( 2 ) 


(3) 


on Bp at z = 0. Here Bp and Bp indicate rigid and flexible boundaries, 
respectively; 3p/3n is the pressure derivative normal to the wall surface; w 
is the displacement of the flexible wall in the z-direction; and a dot indicates 
a time derivative. For the analysis considered in this study, the wall displace- 
ment w is determined independently. The acoustic pressure in the enclosure is 
not assumed to affect the wall. Reference 13 has shown this approach to be valid 
for deep enclosures. 

To solve equation ( 1 ) the pressure can be written 


oo oo 

P = ^ ^ (Cp)ij(z, t) Xjjg(x,y) 

i=0 j=0 


(4) 


iirx jiry 

where X^jg = cos cos are the modes of a cavity with hard walls at 

J a b 

x □ 0,a, and y = 0,b. When the flexible wall motions in terms of these, modes 
are expanded and orthogonality is used, 


OO OO 

-p a w = Gij(t) Xjjo 

i=0 j=0 


In this case, 

Gij(t) 


e i j f a f b •• 

— ) ) PawXijo dx dy 

ab Jq 


(5) 


( 6 ) 
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and 


e ij = 



(i = 0, j = 0) 
(Either i f. 0 or j j? 0) 
(i / 0, j t 0) 


(7) 


Equations (1) to (4) demonstrate that boundary conditions for the pressure 
coefficients (Cp)jj(z,t) are not homogeneous. A direct application of separa- 
tion of variables technique for (Cp)^j(z,t) will not work, and a different 
method needs to be adopted. This solution method can be achieved by transform- 
ing a homogeneous differential equation with nonhomogeneous boundary conditions 
into a problem consisting of a nonhomogeneous differential equation with homoge- 
neous boundary conditions (ref. 14). Using equations (1) to (6) and the 
expression 


( c p)ij = ®i j + Z ( z >t) 


( 8 ) 


where are solutions of the associated homogeneous problem and Z(z,t) is 

the solution on the boundary, gives 


32$ ij 1 •• 


3z 2 


— — $ •? a — A-! -i i — A. -?Z + — Z — 


ij 






92z 

3z 2 


(9) 


and 


3$ 


ij 3Z 

+ — = G-; 4 

3z 3z 3 


( 10 ) 


at z = 0 where 

Ajj = (iiT/a) 2 + (jir/b) 2 


( 11 ) 


3 *ij 

Since ■; are the solutions with homogeneous boundary conditions = 0 

XJ 3z 

at z = 0, equation (10) reduces to 


3Z 


ij 


3z 


= G 


ij 


( 12 ) 


at z = 0. Equation (12) is the boundary condition that must be satisfied by 
function Zjj. Any continuous function which satisfies equation (12) is a suit- 
able function for Zjj (ref. 10). Such a function is a polynomial of the form 

Z ± j = (z - z 2 /2d)Gij (13) 
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The associated homogeneous boundary value problem has a solution 


®ii = XI Bijk(t) cos (14) 

k=0 d 

Then the sound pressure distribution inside the enclosure may be determined by 
combining equations (4), (8), ( 1 3 ) > and (14). The result is 

OO CO CO oo oo 

p(x,y,z,t) = ®i jk( t) ^ijk + X Z-^j(z,t) X jq (15) 

i=0 j=0 k=0 i=0 j=0 

The series solution given in equation (15) is uniformly convergent everywhere 
including the flexible boundary (ref. 11). The equation of modal coefficients 
Bijk can be obtained by substituting equation (15) into equation (9), multiply- 
ing by an orthogonal eigenfunction, and integrating over the volume. The result 
is 

®ijk + jk^ijkBijk + w ijkBijk = (16) 


where the equivalent damping in the enclosure (due to wall absorption and vis- 
cous air damping) is included’ through the modal damping coefficient a ijk- The 
modal frequencies in the enclosure w ijk can be obtained from 


1 /2 

^ijk = c£(iir/a)B + (jiT/b)^ + (kir/d)Bj 



where 

(k = 0) 
(k f 0) 



(17) 


( 18 ) 


(19) 


The solutions for the panel motions w are included in equation (18) through 
the term 


When initial conditions on the coefficients Bijk are assumed, the solu- 
tion to equation (16) can be obtained in the time domain. The linear character 
of the governing equations considered in this study, however, makes it more con- 
venient to solve these equations in the frequency domain. Adopting the Fourier 
integral representation of Bijj c (t) gives 
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1 r * 00 

B ijk( fc ) = iir J 00 6 ~ Wt dW 

and 

Bijij(w) = C BijifCt) e _ i wt dt 

^_oo 


Then the solution to equation (16) in the frequency domain is 


(20a) 


(20b) 


B ijk = 



^ijk 

w2 + 2ia ijk o) ijk o) 


( 21 ) 


where ^ijicC 10 ) are Fourier transforms of F^jj c (t). Since F^jk are functions 
of panel response w, where w is. a Fourier transform of w, the response char- 
acteristics of viscoelastic sandwich panels are determined next. 


RESPONSE OF THE VISCOELASTIC SANDWICH PANELS 

The rectangular sandwich plates considered in this analysis are assumed to 
be flat and simply supported on all four edges. Acting on the top surface of 
the plate is a random noise pressure, and at the bottom face the plate is backed 
by an acoustic cavity as shown in figure 1 . The sandwich plate consists of two 
isotropic elastic plates and an isotropic viscoelastic core. The response anal- 
ysis of the sandwich panel is considered separately for the soft core (compress- 
ible) and for the hard core (incompressible). 


Soft Core 

When the viscoelastic core is very soft, Poisson's ratio of the material is 
nearly zero; consequently, such a material can be approximated by a viscoelastic 
spring. Assuming small-deflection theory, the governing equations of motion are 
(ref. 15) 

V 1 ■■ 

P3h3 + p-|h-| wi + - P3I13W2 + c-|W-| + k3(w-] - W2) = p r (x,y,t) (22) 

1 6 

and 

P3h3 + P2h2jw2 + g P3h3w-j + C2»2 + ^3(^2 - w-j ) = 0 (23) 

where V 1 * = 3V3X 1 * + 2 sVSx 1 * 9y^ + sVSy 1 * and the subscripts 1, 2, and 3 
refer to the top plate, bottom plate, and the core, respectively. The terms 
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3 p 3 h 3 


and 


1 

- P 3 h 3 


represent the apportioned contributions of the mass of the 


viscoelastic core to the displacements w-| and w 2 . The viscoelastic spring 
constant of the core material is k 3 = E 3 (1 + j.3)/h 3 where E 3 is the stiff- 


ness modulus in compression and 8 is the loss factor in the core. In obtain- 
ing equation (23), the effect of cavity pressure on the bottom face plate motion 
has been neglected. 


Solutions to equations (22) and (23) can be represented in terms of the 
simply supported plate modes 

CO CO 

w-|(x,y,t) =2Z2Zc mn (t) 4>mn( x »y) (24) 

m=1 n=1 

GO CO 

w 2 (x,y,t) =^2Z D mnU) <t>nm(x,y) (25) 

m=1 n=1 

where and are the generalized coordinates of the top and bottom 

plates, respectively, and 4> mn = sin (mirx/a) sin (niry/b). Substitution of equa- 
tions (24) and (25) into equations (22) and (23) and use of the orthogonality 
principle give a set of coupled differential equations in and D^. Tak- 

ing the Fourier transform of these equations shows that 

C mn = 4nQmn ( 26a > 


^mn ” ^mnQmn 


(26b) 


Ht 

mn 


7mn 

a 1 (^mn^mn “ °1°2) 


(27a) 


Hb 


mn 


°2 

a l( x mnYmn ~ ° 1 ° 2 ^ 


(27b) 


in which and are the transformed deflection amplitudes of and 

Dmn which correspond to the top and the bottom plates, respectively, are 

the generalized random forces 

Qmn = — f f P r (x,y,to) (j^ dx dy (28) 

ab ^0 
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in which p r are the transformed pressures p r , and 


^mn = + 


Ymn = 


1 / 1 


i — - a) + — ^D 1 [(mir/a) 2 + (rnr/b) 2 ] + k^^j 

+ i — w + — ^ D2^(mir/2) 2 + (mr/b) 2 J + k3^ 


ai = ^\6 p 3 h 3 w + k 3 j 


a 1 

02 - — O' i 
a 2 


a 1 = - P 3 h 3 + Pl h 1 


1 

a 2 = " P3I13 + P2 h 2 


( 29 ) 


(30) 


The natural frequencies of the coupled system can be determined by setting 
c-| = C2 = $ = 0 and using 

^mnYmn ~ Q "\0 2 - 0 (31) 


Equation (31) gives two characteristic values for each set of modal indices 
(m,n). These roots are associated with in-phase flexural and out-of-phase 
dilatational vibration frequencies of the sandwich construction. The dilata- 
tional vibration frequencies are strongly dependent on the core stiffness 
kj = £3/113. For large values of core stiffness k3 these frequencies would 
become very large, and the approach developed in this study would lose practical 
significance. 


The frequency response function for the bottom face plate given in equa- 
tion (27b) can now be combined with equation (21) and subsequently with equa- 
tions (6), (13), and (15) to obtain the pressure distribution inside the enclo- 
sure. Such an expression will be given in the later sections. 


Hard Core 

When the core material in a sandwich construction has a certain degree of 
stiffness in bending and extension, energy is dissipated in the core by the 
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flexural vibrations of the plate. The governing equations of motion for this 
case are taken from reference 2. It is assumed that the normal deflections of 
both face plates are always in phase, the shear strains across the depth of the 
face plate are small, the stresses in the core material parallel to the plate 
surface are negligible, and the in-plane inertia forces are small. The govern- 
ing vibration equations for such a sandwich panel can then be written as 
(ref. 2 ) 


DV 4 w - 

G(h 

+ II 3 ) 

(h + 

h 3 )| 

|, 3 2 w 

2h 3 


E 

3 2 u 

X 

E 

3 2 u 

2G 

1 - V 2 

3x 2 

2(1 

+ V) 

3y 2 

” iih^ 

E 

3^v 



E 

3 2 v 

2G 

1 - V 2 

3x2 

2(1 

+ V) 

3y 2 

hh^ 


' 9u 3v 
(3x + 3y / 


y 3 2 w 1 

' o — ~ = r P r (x,y,t) 

2 3t 2 2 


(32) 


U 


E(3v 2 +1) 32 v Q 3w 

+ (h + ho)— = 0 

V.U. ^ 3 x 


hh 3 2(1 - V 2 )(i _v) 9x 9 y hh 3 


(33) 


where 


E(3v 2 + 1 ) 9 2 u G 3w 

v + — + (h + ho) — = 0 


hl ?3 2(1 - v 2 )(i - v) 9x 3y hh 3 


9 y 


y = 2 ph + p 3 h 3 
G = G 0 ( 1 + iB) 


(34) 

(35) 

(36) 


The displacements u, v, and w correspond to the x-, y-, and z-coordinates, 
respectively, and Gq is the shear modulus in the core. In this formulation it 
is assumed that the thickness h of both face plates is the same. 

For simple support boundary conditions the modes of free vibration can be 
represented by 


u 


v 


w 



^mn 

^mn 

w mn 


rmrx 

cos sin 

a 


IOTTX 

sin cos 

a 


nrrrx 

sin sin 

a 


mry 

niry 

~ 

mry 

~b~ 


(37a) 

(37b) 

(37c) 
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Substituting these equations into equations (32) to (3*0, utilizing the orthog- 
onality principle, and taking the Fourier transformation, give (ref. 2) 


Hnn ~ HjunQjun (38) 

The generalized random forces are defined in equation (28). are the 

Fourier transforms of W^, and 


H 


mn 


1 


Kmn^ 1 + ihmn^ - 


(39) 


^mn 


’’'mn 


'i'mn 



3(1 + t) 2 (1 + Vi + S 2 ) 

f - j n 4 (ri 2 + 1 ) 2 

1 + ■ " 1 1 — 

L (1 + Vi> 2 + 6 2 J 


38 W 1 + t)2 

(1 + 4^) 2 + ^ 2 + 3(1 + t ) 2 (1 + 4 ^ + g 2 ) 

ri 2 + 1 ff2 E / hn \ 2 

(1 _ v) 2 2 G 0\ b / 


T = 


I13 

h 


(40) 

(41) 

(42) 

(43) 


bm 

r 1 = — 
an 


(44) 


The natural frequencies of the sandwich plate can be obtained from (ref. 2) 


f "“ = a 


yb 1 * 


-. 1/2 


1) 2 K, 


mn 


(45) 


The frequency response function for the vibration of a sandwich panel with 
a hard core given in equation (39) can be combined with equation (21) to deter- 
mine noise pressure inside the enclosure in a fashion similar to that described 
earlier for a panel with soft core. 

For the analysis presented in this paper, it is assumed that the input 
pressure spectral density is specified. Thus, the response (panel deflection 
and noise pressure inside the enclosure) needs to be expressed in the form of a 
spectral density. Following the procedure given in reference 16 shows that the 
spectral density of vertical deflections w can be determined from 
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S w (x,y,co) 





H mn H rs4 ) mn ( } ) rsSmnrs 


(46) 


where are given in equations (27a), (27b), and (39), S mnrs are the cross - 

spectral densities of the generalized random forces Qmn> anc * a star indicates 
conjugation. If the input pressure p r is assumed to be a stationary random 
process, using equation (28) and reference 16 gives 

16 P a f» b p a ^b 

S mnrs = S r (C,n,u) (fmntx-i ,y-| ) <l> rs (x2,y2) 6*1 dx 2 dy\ dy 2 (47) 

a 2 b 2 ^0 ^0 ^0 

where S r (£,u,u)) are the cross-spectral densities of the input pressure p r , 
and £ = x 2 - xi , U = y 2 - yi are the spatial pressure separations. For the 
cases where the random process p r has similar spatial correlation properties 
at all delays, the input pressure cross-spectral density can be written as 

S r (£ ,n ,u) = R x (S,w) R y (n,to) S p r (w) (48) 

in which R x and Ry are the narrow band spatial correlation coefficients cor- 
responding to the x- and y-coordinates, respectively, and Sp r is the speci- 
fied input pressure spectral density. 

The analysis developed in this paper for sandwich panel vibrations consid- 
ered separately a soft core and a hard core. A general procedure which would 
not limit the analysis to a particular core stiffness can be developed using 
Lagrange’s equation and the expressions for strain and kinetic energies. An 
attempt to solve such a problem has been presented in references 8 and 9. How- 
ever, as pointed out in reference 7, there remain some significant uncertainties 
in these formulations and solutions. Thus, an in-depth study on the vibration 
characteristics of viscoelastic sandwich panels with no limitations on core 
stiffness would be a useful extension of the work presented in this paper. 


ACOUSTIC-STRUCTURAL MODEL 


The equations developed in previous sections can be combined to construct a 
noise transmission model for a viscoelastic sandwich panel. Applying a Fourier 
transform to equation (15) and using equations (13) and (21) give 


p(x,y,z,co) = 


F ijk x ijk 


i=0 j=0 k=0 a) 2 jk - a) 2 + 2ia ijk w ijk u> 


- z 2 
2d 


1=0 j=0 


GijXijO (49) 


Equations (6), (25), and (26b) demonstrate that 




^mn^mn 


(50) 
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When equations (18), (49), and (50) are combined, the pressure inside the enclo- 
sure caused by the vibration of a viscoelastic sandwich panel with a soft core 
can be written as 


p(x,y,z,w) = 


p a w2 c 


oo oo oo 


oo oo 


ab i=0 j=0 k=0 


ijf"®ijk^ijk + e k Z d( z ^l ^ b ijmn H mnQmn x ijk 

m= 1 n= 1 

(51) 


where 


Z d (z) 


9 ij0 = 
9 ijk = 


H ijk = 



1 

w iik - w2 + 2 i a i jk^ijk^ 


(52) 

(k = 0) (53a) 

(k j. 0) (53b) 

(54) 


■'ljmn 


n b 

x ijo( x >y) 4 >inn( x »y) dx d y 

I 


(55) 


Noise pressure inside the enclosure caused by the vibration of a viscoelastic 
sandwich panel with a hard core can be obtained from equation (51) by replacing 
H b with Hmn where the latter expression is given in equation (39). 


The spectral density of the interior pressure p can be obtained by taking 
the mathematical expectation of equation (51) and then following the procedure 
presented in reference 16. The result is 


S p (x,y,z,w) = 


f p a uj2 c \ 


ab 


i, j,k,r,s,q 


e ij e rsi 


[ 9 ijk H ijk + «k z d( Z )]j 


9 rsq H rsq 


e q z d 


( Z )J b i jmn^mn x i jk^mnlo^ b rslo^lo x rsq^* 


m,n,l,o 


(56) 
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where each index indicates a separate summation and a star denotes conjugation. 
For the cases where acoustic damping in the enclosure is small and the acoustic 
modal frequencies are well separated, the contribution to equation (56) from the 
cross terms + e k Z d )L ± j m Xi j k (e rsq H rsq + e q Z d )«L rslo X rsq is relatively 

small. It can therefore be neglected. When the structural damping is small, 

similar restrictions can be imposed on the cross terms L^j^H^Sjnn^LpgiotHiQ)*. 
However, for viscoelastic sandwich panels, damping is usually large and no such 
simplifying assumption is valid. 


RESULTS AND DISCUSSION 

The numerical results presented in this paper correspond to the acoustic 
cavity and the sandwich panel shown in figure 1 . The cavity beneath the panel 
is assumed to be acoustically sealed. For the calculations made, the top and 
bottom face plates are assumed to be aluminum alloy whereas the core is assumed 
to be a lightweight low modulus viscoelastic material. The physical data used 
are given in table 1 . The input random pressure p r acting on the top face of 
the flexible sandwich panel is taken to be that of truncated Gaussian white 
noise for which the spectral density is 


V = 


0.00058 (N/m 2 ) 2 /Hz 
0 


(0^ f i 1000 Hz) 
(Otherwise) 


(57) 


The spatial pressure distribution is assumed to be uniform over the panel sur- 
face. The spectra given in equation (57) correspond to a 100-dB level at the 
specified frequency range and to about a 130-dB level overall. 

The modal damping in the enclosure was taken as 

a ijk = a 00l( w 001 /w ijk) 2 < 58 > 

where the fundamental modal damping coefficient is assumed otgoi H 0*01 (hard 
wall cavity) or <*001 = 0*10 (cavity with absorbing walls). The acoustic modal 
damping behavior given in equation (58) was observed experimentally in refer- 
ence 14 for a hard wall cavity. The effect of acoustic wall absorption on the 
interior noise in a cavity is strongly influenced by frequency, wall geometry, 
and the type of absorbing material. However, for small acoustic absorption, 
equation (58) is a useful approximation of modal damping for low frequencies. 

The deflection response spectral density at the middle of a sandwich panel 
with a soft core is presented in figure 2 for several values of core loss factor 
8. The peaks in the spectral density correspond to flexural and dilatational 
modes. The fundamental dilatational mode frequency is 401 Hz. Figure 2 shows 
how the low frequency is dominated by flexural modes, whereas at mid-range fre- 
quencies (400 to 600 Hz), dilatational modes are also excited for low values of 
core loss factor 8. Since damping in the face plates and in the core material 
was taken to be constant at all frequencies, the response levels at frequencies 
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above 650 Hz are very low and no distinct peaks are observed. For the case 
where the core damping factor 8 is zero, the core acts as an elastic material 
with no damping and the response at the first dilatational mode is very large. 
The noise transmission characteristics for the same panel are shown in figure 3- 
The noise transmission curve was obtained from 

S p (f) 

NT(f ) = 10 log (59) 

S p r (f) 

where Sp is sound pressure spectral density in the interior calculated from 
equation (56). The input pressure spectral density Sp r is defined in equa- 
tion (57). The significant structural panel modes and the acoustic cavity modes 
are indicated in figure 3. These modes were calculated from equation (31). As 
can be observed from figure 3 , the noise transmission characteristics of an 
elastic panel and a viscoelastic panel with soft core are similar for frequen- 
cies up to about 250 Hz and above 600 Hz. In the intermediate frequency range, 
structural modes corresponding to dilatational vibrations have a significant 
effect on noise transmission. At the frequency of the first dilatational mode 
(401 Hz), a panel with a viscoelastic core can achieve about 30 dB more noise 
reduction. These results indicate that noise transmission into the enclosure is 
dominated by flexural modes for frequencies below 250 Hz, by dilatational and 
acoustic modes in the frequency region 250 to 500 Hz, and by acoustic modes 
above 500 Hz. 

For the viscoelastic sandwich panel with a hard core, damping in the face 
plates was assumed to be negligible in comparison to the damping in the core. 

The thickness of each metallic face plate was taken to be equal to 0.00051 m. 
Viscoelastic cores with two different values of shear modulus Gq were selected 
for the study. For each of these cases, the core loss factor 8 and the core 
thickness I 13 were varied. The displacement response spectral densities at the 
center of the panel are shown in figures 4 and 5 for several values of core 
damping factor 8 . For low values of core loss factor 8 , distinct peaks can be 
observed at the natural frequencies of the panel. The natural panel frequencies 
were calculated from equation (45). Since the core stiffness shown in figure 5 
is much larger than that shown in figure 4, the peaks associated with modal fre- 
quencies are shifted to the right in figure 5. The deflection response is lower 
for the panel with a stiff er core. These calculations were based on core thick- 
ness hg a 0.00635 m. Similar results are presented in figures 6 and 7 for dif- 
ferent core thicknesses. These figures illustrate that response spectral densi- 
ties are significantly affected by change in core thickness and core stiffness. 
Viscoelastic panels with thin but stiff cores respond more than panels with less 
stiff cores. However, when the core thickness is about 0.008 m or larger, panels 
with stiff cores respond less than panels with not very stiff cores. Noise 
transmission characteristics for these panels are given in figures 8 to 13 for 
an enclosure with low (agoi = 0.01) and high (c*001 = 0.1) interior acoustic 
damping. When the acoustic damping in the enclosure is high, the peaks at the 
cavity modal frequencies are suppressed. The two acoustic damping coefficients 
chosen could be representatives of a cavity with hard or absorbing walls, 
respectively. For low acoustic damping and large core loss factor 8 , noise 
inside the enclosure is dominated by the acoustic cavity modes. Approximately 
20 dB more noise reduction can be achieved at the first panel mode by increasing 
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the core loss factor from 0.1 to 1.0. Furthermore, by increasing cavity modal 
damping from 0.01 to 0.1, about 10 dB more noise reduction can be gained at most 
acoustic modal frequencies. Figures 12 and 13 show that a significant amount of 
noise reduction can be achieved by increasing the core thickness. For example, 
by increasing the core thickness 10 times, about 20 and 15 dB more noise reduc- 
tion is obtained at the first panel mode for a core with Gq = 68 900 000 and 
Gq = 2 760 000 N/m^, respectively. In figure 14 the overall sound pressure 
levels in the interior are plotted against the core thickness. The increase in 
total panel mass with the increase in core thickness is also shown in figure 14. 
The mass curve is obtained under the condition that the density of the visco- 
elastic core material is one-fifth the density of the face plates. The overall 
sound pressure was obtained from 

OASPL = 10 log (a2/p Q 2) (go) 

where is the variance of the interior pressure and pq is the reference 

pressure pq = 0.00002 N/m^. These results indicate that the overall noise lev- 
els are about 5 dB higher for a panel with a stiff core (Gq = 68 900 000 N/m^). 
By increasing ' the core thickness 10 times, the overall noise level was reduced 
by about 12 dB. 

A direct comparison between elastic and viscoelastic panel response is 
shown in figure 15. The curves corresponding to the viscoelastic sandwich panel 
were obtained for 8 = 1.0 and hg □ 0.00635 m. The thickness of the elastic 
panel was taken to be equal to 0.00102 m, and the structural damping character- 
istics were represented by 

“ran = a 1 1 1 1 / w mn ) (61) 

where are the modal coefficients and are the natural frequencies of 

the elastic plate. The results shown in figure 15 correspond to i = 0.02. 

The mass of the elastic panel was adjusted to be equivalent to the mass of the 
viscoelastic sandwich panel. A comparison of noise transmission by elastic and 
viscoelastic panels is shown in figure 16. As can be observed from figures 15 
and 16, panel response and noise transmission for an elastic panel is strongly 
dominated by structural vibrations while the noise transmission by a viscoelas- 
tic sandwich panel is dominated by acoustic cavity modes. At low frequencies 
(below 200 Hz) noise reductions of about 50 dB or more can be achieved by visco- 
elastic sandwich panels. However, at frequencies where the acoustic modes are 
dominant, no significant noise reduction is achieved by viscoelastic panels. 

The overall sound pressures in the interior are plotted against core thickness 
for equivalent (mass) elastic and viscoelastic panels. (See fig. 17.) A direct 
comparison of these results indicates that for a viscoelastic panel with a thin 
core (ho = 0.0025 m), the overall noise in the cavity is about 5 dB 
(Gq = 68 900 000 N/m^) and 10 dB (Gq = 2 760 000 N/m*) lower than the overall 
noise for an equivalent elastic panel. For a thick core (hg = 0.025 m) , these 
noise reduction values are 13 and 17 dB, respectively. The higher noise reduc- 
tion values for a thicker core can be attributed to an increase in the sandwich 
construction stiffness and a larger capacity to dissipate vibrational energy. 
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CONCLUDING REMARKS 


An analytical study was conducted to determine noise transmission charac- 
teristics of viscoelastic sandwich panels. The results indicate that noise 
transmission by sandwich panels is strongly dependent on thickness, damping, and 
material properties of the viscoelastic core. 

Sandwich panels with very soft viscoelastic cores transmit noise much like 
elastic panels except in the frequency range where dilatational (out-of-phase) 
modes are excited. About 30 dB more noise reduction can be achieved by visco- 
elastic sandwich panels in this frequency range (300 to 600 Hz). 

The vibration response and noise transmission of sandwich panels with hard 
cores are low when compared with equivalent elastic panels. As much as 50 dB 
more noise reduction can be achieved by viscoelastic panels at some frequencies 
in the low frequency range (below 200 Hz). At frequencies above 200 Hz, acous- 
tic modes dominate interior noise for the acoustic enclosure chosen in this 
study. About 10 dB more overall noise reduction cap be gained by increasing 
the core thickness 10 times (from 0.0025 to 0.025 m) . With increasing core 
stiffness (from a shear modulus of 2 760 000 N/m^ to a shear modulus of 
68 900 000 N/m^), noise reduction decreased by about 5 dB overall. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 11, 1977 
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TABLE 1.- PHYSICAL DATA USED IN THE STUDY 


Enclosure 

a = 0.250 m, b = 0.508 

m, d = 0.762 m 

(fig. 1 ) 

P a = 1.225 kg/m3, c = 

330 m/sec 

Location in the cavity 



where interior noise 

x = 0.102 m, y = 0.152 

m, z = 0.254 m 

was computed 




Soft core 

Hard core 


E-) = E 2 = 72 400 000 000 N/m 2 

E = 72 400 000 000 N/m 2 

Elastic face plates 

ci = c 2 = 27.13 N-sec 2 /m3 

O 

II 

c\ 

0 

II 

6 " 

(aluminum) 

hi b h 2 = 0.00051 m 

h = 0.00051 m 


V-j = v 2 = 0.3 

v = 0.3 


p-l = p 2 = 2770 kg/m3 

p = 2770 kg/m3 


E 3 = 34 500 N/m 2 

G 0 = 2 760 000 N/m 2 



G 0 s 68 900 000 N/m 2 

Viscoelastic 

h 3 = 0.00635 m 

h 3 = Varied 

core 

v 3 = 0 

v 3 = 0.17 


P 3 = O.lp 

P 3 = 0.2p 


O 

t — 

O 

O 

II 

CQ 

6 = 0 . 1 , 0.5, 1.0 

— 
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Figure 1 . - Geometry of enclosure 



(b) Viscoelastic sandwich panel. 


viscoelastic sandwich panel. 
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Figure 3.- Noise transmission by viscoelastic panel with soft core. 

= 0.00635 m; otQOI = 0.01. 





400 

FREQUENCY, Hz 

Figure 5.- Displacement response spectral density for panel with hard core and large shear modulus. 

A = 0.000645; h 3 = 0.00635 m; G 0 = 68 900 000 N/m 2 . 




Figure 7.- Displacement response spectral density for different hard core thicknesses. 
Large shear modulus. 3 = 0.5; A = 0.000645; Gq = 68 900 000 N/m^. 
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Figure 8.- Noise transmission by viscoelastic panel with hard core and small shear modulus. 

ct 001 = 0.01; h 3 = 0.00635 m; G 0 = 2 760 000 N/m 2 . 
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Figure 9.- Noise transmission by viscoelastic panel with hard core and large shear modulus. 

ctQoi = 0.01; h 3 = 0.00635 m; G 0 = 68 900 000 N/m 2 . 
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Figure 10.- Noise transmission by panel with hard core into cavity with large acoustic damping 
Small shear modulus, hg = 0.00635 m; Gq = 2 760 000 N/m^; otQO 1 = 0*1 • 
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Figure 11.- Noise transmission by panel with hard core into cavity with large acoustic damping. 
Large shear modulus, h^ = 0.00635 m; Gq = 68 900 000 N/m^; otQ0 1 = 0*1* 
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Figure 12.- Noise transmission by viscoelastic panels with different hard core thicknesses into cavity 
with large acoustic damping. Small shear modulus, aooi = 0.1; 3 = 0.5; Gq = 2 760 000 N/m^. 
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Figure 13.- Noise transmission by viscoelastic panels with different hard core thicknesses into cavity 
with large acoustic damping. Large shear modulus. Gq = 68 900 000 N/m 2 ; 8 = 0.5; agoi =0.1. 
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Figure 15.- Displacement response spectral density for elastic and viscoelastic panels with hard cores. 

a-|-| = 0.02; 3 = 1.0; A = 0.000645; h 3 n 0.00635 m. 
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Figure 16.- Noise transmission by elastic and viscoelastic panels with hard cores. 

hj = 0.00635 m; 3=1.0. 
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Figure 17.- Overall sound pressure levels for elastic and viscoelastic panels with varying 

hard core thicknesses, a^-j = 0.02; 3 = 1.0. 
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